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System and Method for Statistical Timing Analysis of Digital Circuits 

FIELD OF THE INVENTION 

5 This invention relates to design automation of digital integrated circuits. More 

specifically, it relates to static timing analysis of digital circuits in the presence of delay 
variations. Yet more specifically, the invention relates to propagating arrival and 
required arrival times of a digital circuit in a probabilistic or statistical fashion, while 
considering and preserving correlations between the delays of various circuit elements. 

10 

RELATED APPLICATIONS 

C. Visweswariah, "System and Method for Probabilistic Criticality Prediction of Digital 
Circuits," Docket number YOR92003-0402US1, U.S. Patent Application Number (to be 
1 5 assigned) filed on 09/1 8/03 . 

C. Visweswariah, "System and Method for Incremental Statistical Timing Analysis of 
Digital Circuits," Docket number YOR92003-0401, U.S. Patent Application Number (to 
be assigned) filed on 09/18/03. 

20 

The descriptions set forth in these co-pending applications are hereby incorporated into 
the present application by reference in their entirety. 
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BACKGROUND OF THE INVENTION 



With each succeeding generation of integrated circuit technology, variability is 
proportionately increasing. The sources of such variability include manufacturing 
5 variations, device fatigue, environmental variations and phase-locked loop (PLL) 
variations. In the case of manufacturing variations, the front-end-of-the-line (FEOL) 
which are the layers that define the active transistors show variation in the transistor's 
electrical characteristics. Physical quantities such as the length of the gate, depth of the 
semiconductor junction or thickness of the oxide cannot be perfectly controlled during 

1 0 manufacturing and hence show variations, which lead to variations in the behavior of the 
transistors. As the physical dimensions get smaller in modern technologies, variability is 
proportionately increasing. In addition, the back-end-of-the-line (BEOL), which consists 
of the metal interconnect layers, also exhibits variability. For example, the thickness, 
width and inter-layer dielectric thickness of each metal layer are sources of variability. 

15 These in turn cause the wires to change their delay, and in fact these sources of variability 
can change the delay of gates which are driving them and gates which are driven by them. 

The second main type of variations is due to device fatigue effects such as hot electron 
and negative bias temperature instability (NBTI). After a long period of use in the field, 
20 transistor characteristics change due to these physical phenomena, leading to changes in 
the delay of circuit components. 
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The third main type of variations is due to environmental effects such as temperature and 
power supply voltage. 

The fourth main type of variations is PLL variations which can include PLL jitter and 
5 duty-cycle variability. 

It is to be noted that in addition to the above, there are other sources of variation such as 
model-to-hardware miscorrelation, silicon-on-insulator (SOI) history effects and coupling 
noise. These other types of variation can also be considered during statistical timing 
10 analysis of digital integrated circuits. 

The variation of delays shown by gates and wires in an integrated circuit can be classified 
in many different ways. The variation may be from batch-to-batch during the 
manufacturing, wafer-to-wafer, chip-to-chip or within a single chip. Lens aberration 

1 5 effects during photolithography, for example, can cause variation of the effective length 
of transistors across a reticle field. There can be temperature and power supply voltage 
variations across a chip. The variations can also be classified by the time scales during 
which variability develops. For instance, fatigue effects cause variability over a period of 
years, whereas across the chip temperature or power supply gradients can develop over 

20 seconds or milliseconds, and coupling noise variations can occur in nanoseconds or 

picoseconds. Whichever way they are classified, it is abundantly clear that these sources 
of variation are making integrated circuit analysis and design more difficult and must be 
accurately accounted for during timing analysis. 
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The traditional timing methodology to handle such variability is to conduct multiple static 
timing analyses at different "cases" or "corners" to determine the spread of performance 
of the circuit under these variations. Corners may include, for example, "best case," 
5 "nominal" and "worst case." Unfortunately, the traditional methodology is breaking 
down because the number of independent and significant sources of variation is 
numerous, and too many timing runs would be required. One way to combat this is to 
worst-case or guard-band against some sources of variation, but this causes pessimism in 
the performance prediction. Another way to combat the explosion of timing runs 
1 0 required is to skip the analysis at certain corners, but this is risky since the performance 
of the circuit may be unacceptable at the skipped corners and this may be manifested by 
chips failing on the tester or in the field. Because of these effects, traditional timing 
methodologies are rapidly becoming burdensome, as well as risky and pessimistic at the 
same time. 

15 

A solution to the problems faced by traditional timing methodologies is statistical or 
probabilistic timing analysis. In such an analysis, timing quantities such as delays, arrival 
times and slacks are not treated as single numbers, but rather as probability distributions. 
Thus the full probability distribution of the performance of the circuit under the influence 
20 of variations is predicted by a single timing run. The problems of unnecessary risk, 

excessive timing runs and pessimism are all potentially avoided. Four examples of such 
statistical timing methods in the prior art include Liou et al [J-J. Liou, K-T. Cheng, S. 
Kundu and A. Krstic, "Fast statistical timing analysis by probabilistic event propagation," 
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Proc. Design Automation Conference, June 2001, Las Vegas, NV, pages 661 — 666], 
Scheffer [L. Scheffer, "Explicit computation of performance as a function of process 
variation," Proc. ACM/IEEE workshop on timing issues in the specification and synthesis 
of digital systems, December 2002, Monterey, CA, pages 1 — 8], Gattiker et al 

5 [A. Gattiker, S. Nassif, R. Dinakar and C. Long, "Timing yield estimation from static 
timing analysis," Proc. IEEE International Symposium on Quality Electronic Design 
(ISQED), 2001, pages 437-^42] and Jess et al [J. A. G. Jess, K. Kalafala, S. R. Naidu, 
R. H. J. M. Otten and C. Visweswariah, "Statistical timing for parametric yield prediction 
of digital integrated circuits," Proc. Design Automation Conference, June 2003, Anaheim, 

1 0 CA, pages 932 — 937]. The references cited above are herein incorporated by reference 
in their entirety. 



PROBLEMS WITH THE PRIOR ART 



1 5 There are several reasons why prior-art methods are not suitable in practice. The main 
shortcoming of prior-art methods is their failure to correctly account for correlations. 
The delays of gates and wires on an integrated circuit exhibit strong correlation. 
Consider a simple example to demonstrate the importance of correlations. A chip that 
has 50,000 latches has 50,000 setup timing tests and 50,000 hold timing tests. Assume 

20 that the probability of each of these tests being successfully met is 99.99%. If the 
100,000 probabilities are perfectly correlated, then the overall probability of the chip 
working correctly is 99.99%. That is, if one test passes, they will all pass. However, if 
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the probabilities are independent, then the probability of making a working chip is 0.9999 
raised to the 100,000 th power, which is a paltry yield of 0.005%! 

There are many sources of delay correlation, and some examples are described below. 
5 Two paths may share some of the same gates; in this case, the delays of the two paths are 
correlated. If a particular chip has stronger P-type transistors due to some manufacturing 
variations, chances are that every single P-type transistor on that chip will be stronger. If 
the 3 level of metal is a little thicker, it is likely to be thicker across the entire chip. A 
launching path (path that gets the data signal to a latch) and a capturing path (path that 

10 gets the corresponding clock signal to the same latch) may exhibit some commonality 

and therefore have correlated delays. The commonality could include sharing some gates 
along the path, sharing metal layers, sharing a power supply voltage island, sharing gates 
of the same type, etc. Delays of gates may also be correlated because of their physical 
position on the surface of a chip. For example, two gates that are close to each other are 

1 5 unlikely to see significant relative channel length variation and their delays are therefore 
likely to be tightly correlated. 

All of these sources of variation must be accurately taken into account, or else the results 
will not be meaningful. In Liou et al, the delay of each individual gate is considered to be 
20 completely independent of any other gate, rendering the analysis unusable in practice. In 
Gattiker et al, gate delay correlations are considered. Critical paths are enumerated, but 
when the effects of these paths are combined, the delay of each path is considered to be 
independent of any other, thus ignoring several important sources of correlation. 
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Among the prior art methods, one that takes into account correlation is Jess et al. This 
method is a path-based method. The delay or slack of each path is collected and 
represented as a first-order model of the sources of variation. Then the slack of the 
5 overall circuit or chip is computed by combining these path slacks in a correlated 
probabilistic fashion. Unfortunately, there are an exponential number of paths in any 
circuit. It is not realistic to list and analyze all the paths in the circuit. Jess et al suggest 
that the top N critical paths be considered, but of course there is no guarantee that the 

(N + l) f / path (or any path other than the first N paths) will not be critical at some point 

10 or corner in the process space or space of variations. All path-based methods have the 
fundamental limitation that the number of paths is too large and some heuristic must be 
used to limit the number of paths submitted for detailed analysis. 

It is to be noted that although there are many significant sources of correlation in the 
1 5 delay variability of integrated circuits, there are some completely random sources of 

variation as well. For example, the oxide thickness in transistors in a modern technology 
is only a few atoms thick and for various reasons it is possible for transistors to have one 
more or one less layer of atoms, leading to variations that are quite random from 
transistor-to-transistor. While the prior-art method of Liou et al can handle such random 
20 variations, other prior-art methods such as those of Gattiker et al and Jess et al cannot. 

Further, all of the above prior-art methods have other problems. An important output 
from a static timing program is diagnostics that can enable a designer, or automated 
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synthesis or optimization program to improve the circuit. Here again, prior-art methods 
are lacking. For example, if a critical path has excessive sensitivity to a particular source 
of variation, then the circuit is not robust in the face of variations. If such diagnostics 
were available, the designer or synthesis program could take various measures to reduce 
5 the excessive sensitivity. So it is important for the output of a statistical timing program 
to be parameterized by the sources of variation. Instead of reporting that the arrival time 
at a node of the circuit has an unacceptably late probability distribution, it would be 
useful to report the sensitivity of the late arrival time to the various sources of variation in 
order to enable suitable remedies. Prior-art methods produce statistical timing results as 
1 0 probability distributions. Unfortunately, this does not help a human designer or 

automated optimization program in improving the performance or robustness of the 
circuit. 

For these and other reasons, statistical timing methods that have been proposed in the 
1 5 literature are not used in industrial practice. 

ASPECTS OF THE INVENTION 

An aspect of this invention is an improved system and method for statistical or 
20 probabilistic static timing of a digital circuit. 
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Another aspect of this invention is a method for statistical timing that has linear 
complexity in the size of the timing graph, and linear complexity in the number of 
sources of variation. 

5 Another aspect of this invention is a method for statistical timing that takes into account 
correlations between delays of individual gates and wires, and correlations between 
delays of paths of the circuit. 

Another aspect of this invention is a method for statistical timing that allows delay 
1 0 models that contain a deterministic part, a correlated random part and an independent 
random part. 

Another aspect of this invention is a timing result that expresses the arrival time, required 
arrival time, slew and slack of each node of the timing graph of the digital circuit as a 
1 5 probability distribution. 

Another aspect of this invention is a timing result that expresses the arrival time, required 
arrival time, slew and slack of each node of the timing graph of the digital circuit in a 
form that is parameterized by the sources of variation. 

20 
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SUMMARY OF THE INVENTION 



The present invention is a system and method for statistical or probabilistic static timing 
analysis of digital circuits, taking into account statistical delay variations. The delay of 
5 each gate or wire is assumed to consist of a nominal portion, a correlated random portion 
that is parameterized by each of the sources of variation and an independent random 
portion. Arrival times and required arrival times are propagated as parameterized random 
variables while taking correlations into account. Both early mode and late mode timing 
are included; both combinational and sequential circuits are handled; static CMOS as 
1 0 well as dynamic logic families are accommodated. The timing analysis complexity is 
linear in the size of the graph and the number of sources of variation. The result is a 
timing report in which all timing quantities such as arrival times and slacks are reported 
as probability distributions in a parameterized form. 

15 BRIEF DESCRIPTION OF THE FIGURES 

The foregoing and other objects, aspects, and advantages will be better understood from 
the following non-limiting detailed description of preferred embodiments of the invention 
with reference to the drawings that include the following: 

20 

Figure 1 is a block diagram of one preferred embodiment of the invention depicting 
statistical or probabilistic static timing of a digital circuit. 
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Figure 2 is a flow chart of the preferred method of conducting statistical or probabilistic 
static timing of a digital circuit. 

Figure 3 shows two edges and three nodes of a timing graph to illustrate the basic 
5 forward propagation operations involved in statistical timing analysis. 

Figure 4 shows two edges and three nodes of a timing graph to illustrate the basic 
backward propagation operations involved in statistical timing analysis. 

1 0 Figure 5 is a block diagram of one preferred embodiment of an output report. 

DETAILED DESCRIPTION OF THE INVENTION 

An inventive statistical or probabilistic static timing flow is shown in Figure 1 . The first 
1 5 input is the netlist representing the structure of the circuit to be analyzed, shown in box 
100. The second input is a set of timing assertions, box 1 10. These typically include 
arrival times at the primary inputs, required arrival times at the primary outputs, 
information about the phases of the clock, and details of external loads that are driven by 
the primary outputs. The assertions can be in the form of deterministic numbers or 
20 independent probability distributions or correlated probability distributions. The third 
input is a set of parameterized delay models, box 120. These allow the timer to 
determine the delay of a gate or wire as a function not only of traditional delay-model 
variables (like input slew or rise/fall time, and output load) but also as a function of the 
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sources of variation. For example, a first-order linear model may be employed, like the 
one shown below: 



where the delay consists of a deterministic (constant) portion a 0 , a correlated (or global) 



5 portion Va,Axy and an independent (or local) portion a n+} AR . The number of sources of 



variation is n , and a j , i = 1, A , n are the sensitivities of the delay to the sources of 
variation x i , i = 1, A , n and a nJrX is the sensitivity to an independent random source of 
variation R . The notation Ax ; denotes the deviation of x i from its mean or nominal 

value, and AR denotes the deviation of R from its mean or nominal value. It is to be 
1 0 understood that the delay models can be stored in a pre-characterization step or calculated 
on the fly as required. The format in which they are stored could include analytical delay 
equations or table models. The next input is information about the statistics of the 
sources of variation, box 130. This input typically has a list of the sources of variation 
with a mean value and standard deviation for each source of variation. Any correlations 
1 5 between the sources of variation are specified here. 

Circuit netlists 100, assertions 110, parameterized delay models 120, statistics of the 
sources of variation 130, and their functional equivalents are well known in the prior art. 

20 The novel probabilistic or statistical static timing program, box 140, accepts all of these 
inputs and produces a novel statistical timing report, box 150. This report typically 



n 




n 
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includes arrival times, required arrival times, slacks and slews at each node of the circuit. 
These timing quantities are not single numbers, but rather they are probability 
distributions. The information in the timing report can take many forms, including one or 
more of: mean value and variance for each timing quantity, a parameterized 
5 representation of the distribution of each timing quantity, a graphical representation of 
the distribution of each timing quantity, and a correlation report between these various 
timing quantities. Various automatic audits such as checking for excessive sensitivities 
can be built into the timing report. The excessive sensitivities are of interest since they 
must be reduced in order to improve the robustness of the circuit. It is to be understood 
10 that the circuit being timed can be very large and can consist of millions of gates and 
wires. All the information mentioned above could be extremely voluminous, and so it is 
ordinary practice to provide options to selectively output the required information, or 
even calculate or graphically display all this information on-demand. 

* 

15 The details of the novel statistical timer of box 140 are shown in flow 200 of Figure 2. 
The first step is to read the netlist that contains details of the topology of the circuit to be 
timed. Since this netlist is often in hierarchical form, it is flattened (i.e., the number of 
levels in the hierarchy is reduced). The assertions (each of which can be either 
deterministic or probabilistic) are read and so are the parameterized delay models. 

20 Information about the sources of variation such as the mean and standard deviation 

values and any correlation between the sources of variation are stored in memory in box 
210. 
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The next major step, shown in box 220, is the construction of the timing graph, which is a 
step that is familiar in all static timing programs. In the graph, each node represents a 
node or signal of the circuit and each arc or edge represents a delay in the circuit incurred 
when a logical transition (from low to high or high to low) is transmitted through a circuit 
5 component such as a gate or wire. All possible valid logical transitions of the circuit are 
therefore captured in this graph. Arrival times are typically stored on the nodes of the 
graph and delays of individual gates and wires on the edges of the graph. An arrival time 
in late mode is the earliest time at which the corresponding signal is guaranteed to be 
stable at its correct logical value (having gone through any of the possible paths of the 

10 electrical circuit), and an arrival time in early mode is the earliest time at which the 
corresponding signal can change from its previous cycle stable logical value (i.e., the 
output can not change earlier than the early mode arrival time). Sequential elements and 
dynamic circuits in the graph are represented by a special kind of edge called a test 
segment, which is an indication to the timing program that a timing test must be 

1 5 performed between two nodes of the graph to ensure correct timing operation of the 

circuit. Building of such a graph for both gate-level netlists and transistor-level netlists is 
known in the prior art. 

The next step is an inventive correlated statistical forward propagation of arrival times 
20 through the timing graph, box 230. Next, in box 240, an inventive correlated statistical 
backward propagation of required arrival times is performed. Finally, in box 250, 
inventive timing reports are produced. The details of the calculations involved in boxes 
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230 and 240 for both early mode and late mode analysis are explained in the following 
paragraphs. 

The four basic operations in static timing analysis are "plus," "minus," "max" and "min." 
5 These four basic operations must be replaced by probabilistic equivalents for successful 
statistical timing analysis. It is crucial that the probabilistic equivalents must correctly 
consider and propagate correlations. Consider the simple situation shown in Figure 3 
where two directed edges of the timing graph (340 and 350) from nodes a (310) and b 
(320) meet at a common timing node c (330). The edge labeled d ac represents the delay 

1 0 from node a to node c, and d bc represents the delay from node b to node c. We assume 
that the early and late arrival times at a and b are known to be of the form 



n 




i=\ 



n 




1=1 



n 



AT e a arly = al arly + £«r* A*/ + aZ fy &R e a arly , and 




and the delays are known to be of the form 



n 



late 
S > 



n 




n 



n 



15 
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At node c, it is required to determine AT™ = max[^77 e + d™\ATf e + <f ] and 
42*-* = min ^j^ + d e a f ,AT b early + 1 > which ^ the earliest time at which c is 

guaranteed to be stable in the present clock cycle and the earliest time at which c can 
change from its previous cycle stable logical value, respectively. While the late mode 
calculations are illustrated below, the early mode calculations are performed in a similar 
fashion. Thus, 



ATl ate = max 



n n 



late . V 1 Aate a , late a rtlate . late . Aate a „ . late a relate I 



1=1 7=1 

n n 



= max 



i /ate 

(C +g? e ) + t(°? e +gf' e )Ax,. -W^Atf" f, 

n 



(C +C)+X(^' e H**)Ax, +/£'[A#J 



/=! 



If the sources of variation (the global x and the local R variables) are Gaussian or normal 
distributions, the two quantities above whose maximum we seek to find are also Gaussian. 



1 0 Hence we can write 

AT?* = max[fa = N{^a x )\{p 2 = N(m 2 ,<t 2 )}], 

where p x and p 2 are Gaussian random variables with means and variances assumed to 

be as shown above. These two random variables are not independent, i.e., the two 
random variables are correlated. However, the sources of variation, x, may be correlated 
15 or independent. Each of these cases is described below. 



Assuming that the sources of variation are independent of each other and of zero mean 
and unit variance, we can write down the covariance matrix as 
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cov(/; 1 ,/7 2 ) = 



2,1*/ j + Wn + ij 2/ a / X*i +*/ ) 



+gt' e )(b'r+hf" e ) jrfe- +/^} 2 +fc} 2 



1=1 



1=1 



but if the sources of variation are not independent, the covariance matrix is 



(a,+g,) (a 2 +£ 2 ) A 
fo+A,) (b 2 +h 2 ) A 



(«„+£«) a nM g„ +1 0 0 



(«i + £.) 
M 

(«„ + £„) 
a „ + i 

0 
0 



(Z> 2 + A 2 ) 
M 

0 
0 



where Fis the covariance of the sources of variation supplied by the user and the 
superscript late has been omitted in the latter equation for brevity. (Note that in the 
independent case above, Fwas a diagonal matrix.) In either case, comparing to 



cov(p,,/7 2 ) = 



o x o 2 p 



o x c 2 p a 



2 J 



<r, , a 2 and the correlation coefficient p can be deduced. In the case when the sources of 

variation are independent, this can be achieved by collecting the sum of the squares and 
1 0 products of the corresponding coefficients of the two Gaussians whose maximum we 

seek to determine. It is to be noted that the random components of arrival time and delay 
contribute to the diagonal terms of the 2x2 covariance matrix thus increasing the variance 
of the two quantities whose maximum is being determined, but do not contribute to the 



15 



off-diagonal terms, thereby reducing the correlation coefficient. 



The next step is to compute max[p, ,p 2 ], which can be done by prior art methods as 



taught in [J. A. G. Jess and C. Visweswariah, "System and Method For Statistical 
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Modeling And Statistical Timing Analysis Of Integrated Circuits," U.S. Patent 
Application number 0/184,329 filed with the U. S. Patent Office on June 27, 2002] and 
Jess et al [J. A. G. Jess, K. Kalafala, S. R. Naidu, R. H. J. M. Otten and C. Visweswariah, 
"Statistical timing for parametric yield prediction of digital integrated circuits," Proc. 
5 Design Automation Conference, June 2003, Anaheim, CA, pages 932 — 937]. The 

references cited above are herein incorporated by reference in their entirety. Briefly, the 
preferred way this is performed is to sweep a variable 7] through a suitable range of 

values (for example, from max^ - 3<y x ,// 2 - 3cr 2 ) to max(//, + 3a \ ,ju 2 + 3a 2 ) ). At 

each value of t] , the probability of p x = 7] is multiplied with the conditional probability 

10 of p 2 < t] to obtain the probability that p x dominates p 2 . Likewise, at that same value of 

f] , the probability of p 2 =rj is multiplied with the conditional probability of p x < Tj to 

obtain the probability that p 2 dominates p x . This exercise results in the probability 

distribution of AT l c ate . The individual probabilities of p x dominating p 2 and p 2 

dominating p x are collected over all values of tj ; these are called "binding probabilities" 
15 in the references cited above, whereas in this document they are called "arrival tightness 
probabilities" (ATP) since these are the probabilities that the arrival time at c is 
determined by the upper and lower edges of Figure 3, respectively. Let us denote these 
arrival tightness probabilities as T and (1 - T) , respectively. It is to be noted that the 

method in Jess et al uses a similar technique, but does not consider a purely random 
20 component of the delay, and also applies the method on a path basis, whereas here it is 
being applied at each node of the timing graph. 
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The last step before moving further along in the timing graph is to re-express the late 
arrival time at c in a form that allows it to be propagated downstream in the timing graph. 
The way it is re-expressed is as follows. Since the arrival time at c is determined by the 
upper edge of Figure 3 T % of the time, the dependence of the arrival time at c on the 
5 global sources of variation are derived from the upper edge, but with a weight factor of T . 
Likewise, the dependence of the arrival time at c on the global sources of variation is 
derived from the lower edge, but with a weight factor of (1 - T) . Expressed 

mathematically, 



1 0 All that remains is to determine c* e and c** , which are easily done from the probability 

distribution of AT**' . The mean of the distribution is assigned to c l " e , and the random 

portion c 1 "* is assigned a value such that the variance of the distribution matches the 
variance of the re-expressed arrival time. 

1 5 Now that the arrival time at c has been expressed in the "standard form," it can be 
propagated downstream. It is to be understood that early mode calculations simply 
require computation of the minimum of two distributions instead of the maximum. The 
conditional probabilities are computed in a slightly different way and the rest of the 

procedure is identical. For example, during the sweeping of 7] , the probability of p^-t] 
20 is multiplied by the conditional probability of p 2 >V to obtain the probability that p x is 
smaller than p 2 . It is to be understood that one of ordinary skill in the art can extend the 



n 




T(a'°' e + gf* ) + (1 - T){bf' e + hf e ), i = 1,2, A n. 
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above methods to accommodate various functional forms of the delay equations and 
various types of probability distributions. All that is required is to compute the tightness 
probabilities T and (1 - T) , and to re-express the arrival time at c in the "standard form" 

for downstream propagation. 

5 

When more than two edges converge at a timing node, the maximization or minimization 
is done two edges at a time. The arrival tightness probabilities are stored as the 
computation unfolds, and then the final tightness probabilities are assigned in a post- 
processing step. For example, suppose there are 3 edges that converge at a timing point. 

10 Suppose the tightness probabilities of the first two are determined to be 60% and 40%, 
respectively. Now the maximum of these two distributions is in turn max'ed with the 
path delay represented by the 3 rd edge. Suppose the tightness probabilities that results 
from this computation are 80% for the first pair of edges, and 20% for the 3 rd edge. Then 
the final arrival tightness probabilities for the 3 edges are respectively 48% (0.6x0.8), 

15 32% (0.4x0.8) and 20%. 

It is to be understood that early mode arrival times are propagated forward through the 
timing graph in an analogous manner, the only difference being that early mode delays 
are considered and the "min" operation is used instead of the "max" operation. 

20 

In this manner, both early and late arrival times are propagated forward through the 
levelized timing graph until the arrival times of all end points have been computed. End 
points are either primary outputs or test segments. At primary outputs, the probabilistic 
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arrival time is subtracted from the asserted required arrival time in late mode, and vice 
versa in early mode, to determine the late and early slack of the primary output, 
respectively. A zero slack means that the timing requirement was exactly meant. A 
positive slack means that there is some timing margin over and above meeting 
5 requirements. A negative slack means that the circuit will not perform correctly. 

Similarly, arrival times are compared to each other at test segments to determine a slack. 
For example, at a latch, the late data arrival time is typically added to a setup (or guard) 
time of the latch and then compared to the early clock arrival time to make sure data is 
10 correctly latched. For different types of latches and different types of dynamic circuits, 
these comparisons will take different forms, but they are all straightforward extensions of 
the same concept. 

After the forward propagation is complete, the next step is correlated statistical backward 
1 5 propagation of required arrival times (box 240 of Figure 2). Starting at each end point, 
required arrival times are propagated backwards in a levelized fashion, as in traditional 
static timing analysis. However, in the inventive method, these required arrival times are 
statistical in nature. The operations involved in box 240 of Figure 2 are conducted 
similarly to box 230 of Figure 2, but with a few important differences. The first is that 
20 during backward traversal, delays are subtracted from required arrival times rather than 
being added to arrival times. The second difference is that during backward traversal, 
late mode analysis calls for the "min" operation whereas early mode analysis calls for the 
"max" operation. The third difference is that the required arrival tightness probability of 



YOR920030403US1 



an edge of the timing graph is defined as the probability that the required arrival time of 
the source node of the edge is determined by that edge. These tightness probabilities are 
determined in much the same way as arrival tightness probabilities during the forward 
propagation of box 220, and applied in much the same way to express required arrival 
5 times in the standard manner before propagating further. Thus the backward propagation 
is conducted quite similarly to the forward propagation, but with these few important 
differences. 



Backward propagation is illustrated in reference to Figure 4. The basic goal of backward 
10 propagation is to determine the required arrival time of every node of the timing graph. 
Suppose we are trying to compute the required arrival time at node a (box 410) of Figure 
4, which has two fanout edges (box 440 and box 450) to nodes b (box 420) and c (box 
430), respectively. Then the required arrival time at node a in late mode can be written 
as 



15 RAT l a ate = min[{RAT{ ate -d l £},{RAT l c ate -d 1 ™}] 

and the required arrival time at node a in early mode can be written as 

RA T a early = max[ {RA T b early - d e a f y } , {RA T c earIy - d e a a c rty } ] . 

The actual process of subtracting the delay of the edge and then taking the max or min is 
similar to the operations performed during forward propagation. Likewise, required 
20 arrival tightness probabilities are determined in an exactly analogous manner to the 
corresponding calculations during forward propagation. 
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Once the forward and backward propagation are complete, the late and early arrival times, 
required arrival times and slacks are now available everywhere in the timing graph, and 
these can be reported to the user in many different forms (box 250 of Figure 2). These 
can include graphical plots, textual tables, queries on-demand, audits of failing timing 
5 tests, audits of excessive sensitivity at all timing points or all end points, and so on. 

Timing reports are illustrated in box 510 Figure 5. Timing reports are communicated to 
the user either by means of a programming interface, or by means of a hard disk file or 
files. A timing report typically consists of some circuit information (box 520) and the 

10 corresponding statistical timing information (box 530). Circuit information can include a 
list of gates, components and wires; or a list of paths; or a list of nodes; or a list of 
sequential elements; or a list of end points (primary outputs and timing tests); or a list of 
clock phases. These items can be sorted and filtered in various ways to make the report 
intuitive and productive to the reader of the report. The corresponding statistical timing 

1 5 information, in the case of a node, could include one or more of the node's statistical 
arrival time, statistical required arrival time, statistical slew or statistical slack. For a 
timing test or primary output, the corresponding timing information could include the 
probability that the timing test is met, or the primary output meets its required arrival 
time, respectively. For a path, the corresponding timing information could include the 

20 statistical path slack and statistical arrival time, required arrival time, slew and slack of its 
end point. Further, each statistical timing quantity in the report can be represented in 
various forms, including a mean value and standard deviation; a mean value, independent 
random part and a correlated part; a graphical display of the distribution of the timing 



YOR920030403US1 



23 



quantity; or sensitivities to individual global sources of variation. Further, given any two 
statistical timing quantities, the report could include the correlation coefficient of the two 
quantities, the covariance matrix of the two quantities, and the probability that one is 
larger or smaller than the other. It is to be understood that each of the timing quantities in 
5 the above description can be one of an early-mode or late-mode timing quantity; one of a 
rising or falling timing quantity; and a timing quantity that is specific to a particular 
phase of a particular clock. It is to be further understood that once the statistical timing 
analysis is completed, these results can be reported in a variety of useful ways. 

10 It is to be understood that the detailed description of this invention was explained in the 
context of a simple snippet of a timing graph. One of ordinary skill in the art will be able 
to extend these concepts to accommodate separate rising and falling arrival times and 
delays; sequential circuits; circuits with transparent latches; extensions to handle slew 
(rise/fall time) propagation and effects; probabilistic delay models that are purely random; 

1 5 probabilistic delay models that are correlated; sources of variation that are random; 

sources of variation that are correlated; circuits with probabilistic guard times on the latch 
setup and hold tests; circuits with probabilistic guard times on dynamic circuit timing 
tests; circuits with multiple clock phases; cases in which clock cycle periods and timing 
assertions are probabilistic; transistor level netlists, and situations where the 

20 parameterized delay models are calculated on the fly by circuit simulators. 

Given this disclosure it is apparent to one skilled in the art that the inputs received by the 
probabilistic or statistical static timing process (box 140 of Figure 1) can be any input 
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generally known to computer systems, including but not limited to: keyboard or mouse 
entries, disk, tape, CD-ROM, network connection, fiber optic connection, radio frequency 
link, infra red link, etc. Further the outputs including the parameterized timing report 
(box 150 of Figure 1) can take the form of any known computer output. These outputs 
5 include but are not limited to: printed output from a printer, images on a graphical user 

* 

interface (GUI) or CRT, content on storage media (e.g., memory, CD-ROM, disk, 
diskette), files, information transmitted over a network (fiber optic, telephone, cable, 
radio frequency, infrared, etc.). 
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